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Evidence for particle stability of aaH been suggested by the BNL-AGS E906 experiment. We 
report on Faddeev-Yakubovsky calculations for the four-body AApn system using AA'^ interactions 
which reproduce the observed binding energy of aH(§^) within a Faddeev calculation for the Apn 
subsystem. No aaH bound state is found over a wide range of AA interaction strengths, although 
the Faddeev equations for a three-body AAd model of aaH admit a 1"*" bound state for as weak a 
AA interaction strength as required to reproduce _BAA(AAHe)- 
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I. INTRODUCTION AND INPUT 

Information on hyperon-hyperon interactions is not 
readily available from experiments in free space. It is 
almost exclusively limited to the study of strangeness 
5"=— 2 hypernuclear systems, only handful of which have 
been identified to date. This information is crucial for ex- 
trapolating into multi-strange hadronic matter, for both 
finite systems and in bulk (Ref. and references cited 
therein). Until recently only three candidates, identified 
in emulsion experiments |^, Q , existed for AA hypernu- 
clei. The AA binding energies deduced from these events 
indicated that the AA interaction is strongly attractive in 
the ^5*0 channel Q, in fact considerably stronger than the 
AA^ interaction deduced from single-A hypernuclei, and 
this seemed at odds with the natural expectation borne 
out in one-boson-exchange (OBE) models using flavor 
SU(3) symmetry or within the naive quark model. For 
example, the recent Nijmegen soft-core (NSC97) model 
I, yields 



Vaa ^ Van < Vnn 



(1) 



for the strength V of these essentially attractive interac- 
tions. It is gratifying then that the recent unambiguous 
identification of a ^He in the KEK hybrid-emulsion ex- 



periment E373 yielding binding energy substantially 
lower than that deduced from the older dubious event 
|B|, is consistent with a scattering length aAA~— 0.5 fm 
indicating a considerably weaker A A interaction than 
that specified by oaat'^— 2 fm Q for the AA^ interaction. 
With such a relatively weak AA interaction, and since the 
three-body system AAA^ is unbound (comparing it with 
the unbound Ann system ) , the question of whether 
or not the onset of binding in the S=—2 hadronic sector 
occurs at A— 4 becomes highly topical. 

The Brookhaven AGS experiment E906, studying 
capture following the {K~,K~^) reaction on ^Be, has re- 
cently given evidence for excess pious that defied known 
single-A hypernuclear weak decays and were conjectured 
as due to the formation of ^^H (/=0, J''=l+) A sub- 



sequent theoretical study Q of the weak-decay modes 
available to aa^ does not support this conjecture and 
in our opinion the question of whether or not aa-H is 
particle stable remains experimentally open. If it is con- 
firmed in a future extension of E906 or of a related exper- 
iment, then this four-body system AApn would play as a 
fundamental role for studying theoretically the hyperon- 
hyperon forces as the bound state of the three-body 
system Apn has played for st udy ing theoretically the 
hyperon-nucleon forces (Ref. and references cited 

therein) . Our aim in this Letter is to search theoretically 
for a possible aa-H bound state by solving the appropriate 
Faddeev-Yakubovsky equations for the four-body system 
AApn, particularly for AA interactions which reproduce 
the recently deduced binding energy of aa-Hg @ . This is 
the first ever systematic Faddeev-Yakubovsky calculation 
done for the A=A, S=—2 problem. It has the virtue of 
taking into account properly all the rearrangement chan- 
nels (or equivalently clusters) which the AApn system 
may be split into. We note that a bound state does 
not necessarily imply, for attractive A A interactions, that 
^^H is particle stable. 

The AA^ and AA interaction potentials used as input 
were of a three-range Gaussian s-wave form similar to 
that used by Hiyama et al. jlj, 15 : 



y(2S+i)(^)^^„P+i)exp( 



(2) 



The values of the range parameters (3i and of the singlet- 
and triplet- strength parameters vl are listed in Ta- 
ble |. The AA interaction, respecting the Pauli principle, 
is limited to the singlet s-wave channel. The short-range 
term {i—i) provides for a strong soft-core repulsion and 
the long-range term {i—V) for attraction. The parameter 
7, which controls the strength of the mid-range attractive 
term (i=2), was chosen such that the potential (H) repro- 
duces the scattering length and the effective range for a 
given model as close as possible. Its appropriate values 
for AiV are listed in Table |^ for two versions of model 
NSC97 Q considered realistic ones. For AA we listed 
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TABLE L Range (/3) and strength (v) parameters of the three- 
range Gaussian potential (Gl). 



i 


ft (fm) 


' (MeV) 


(MeV) 


1 


1.342 


-21.49 


-21.39 


2 


0.777 


-379.1 X 7(1' 


-379.1 X 7'^) 


3 


0.350 


9324 


11359 



TABLE II: Values of the parameter 7(23+1) appropriate for 
simulating the Ap potentials of model NSC97 and for a AA 
potential reproducing i3AA(AAHe). The resulting scattering 
lengths a (in fm) are also listed. 



Model 


^(1) 


1 

a 


^(^) 




AN: NSC97e 


1.0133 


-2.10 


1.0629 


-1.84 


AA^: NSC97f 


1.0581 


-2.50 


1.0499 


-1.75 


AA: ^iHe i 


0.6598 


-0.77 







the value which was determined in Ref. Q to reproduce 
the recently reported B\\{j^^^llc) [Q. Also listed are the 
values of the scattering lengths for these AN and AA 
model interactions which obviously satisfy Eq. (||). For 
the pn triplet interaction we multiplied the AiV potential 
(^ by a factor a = 2.0685, using jpn = 1.0498, in order 
to reproduce the NN low-energy scattering parameters 
in this channel plus the binding energy of the deuteron. 
We usedjfor comparison, also the Malfliet-Tjon potential 
MT-III . Our results are insensitive to which form is 
used. 



II. FADDEEV-YAKUBOVSKY CALCULATIONS 



TABLE III: Ba(aH(5^)) and Ad low-energy doublet scatter- 
ing parameters {^Ba in MeV; ^a, in fm) calculated for the 



7 = Apn system. 



4 Dscatt 

-Da 



(in MeV) for ^H(|^) was obtained 



using the effective-range expansion in the quartet Ad channel. 



Model 


''a 






4 Dscatt 

-Da 


NSC97e 


20.7 


2.61 


0.069 


0.015 


NSC97f 


13.1 


2.46 


0.193 


0.003 


NSC97f' 


13.1 


2.46 


0.193 


-0.003 


EFT 


16.8t^;l 


2.3±0.3 


0.13±0.05 




exp. 






0.13±0.05 





the four-body AApn calculation described below. The 



(unobserved and probably unbound) excited state of 
comes out very weakly bound in our Faddeev calcula- 
tion in both versions e and / of model NSC97. In order 
to check the sensitivity of the four-body calculation to 



the location of \ 



H(| ) we will give below results also 

for model NSC97f', where /' coincides with / for the 

channel but slight ly d eparts from it for the ^ channel 
as shown in Table 

Focusing on the AApn Faddeev- Yakubovsky calcula- 
tion we note that for two identical hyperons and two es- 
sentially identical nucleons (upon introducing isospin), as 
appropriate to the /=0, J'^=l+ ground state of ^a-H, the 
18 Faddeev- Yakubovsky components which satisfy cou- 
pled equations reduce to seven independent components, 
in close analogy to the Faddeev- Yakubovsky equations 
discussed in our recent work for the AAaa model of 
^^^Be. Six rearrangement channels are involved in our 
s-wave calculation for aaH: 



We solved the differential Faddeev equations under 
the s-wave approximation jl^ for the I—O, J^^^^ t 
ground-state doublet levels of viewed as a Apn sys- 
tem. Similar calculations for three-body systems are dis- 
cussed in Ref. Some of our results are displayed in 
ground state is bound and the calcu- 



Table im The i 



lated binding energies of the A hyperon (-Ba) are in rough 
agreement with that observed. For model NSC97f, for ex- 
ample, our calculated Ba = 0.19 MeV agrees with that of 
the recent Hiyama et al. [ p8| where no s-wave approxima- 
tion was invoked. The impact of the higher partial waves 
for was estimated by Cobis et al. [|l9j to be of order 
0.02 MeV, well within the error of the measured binding 
energy. Our values satisfy the effective-range expan- 
sion in terms of Ad low-energy parameters which are close 
to those derived using Effective Field Theory methods 
pO| . The convergence of the Faddeev calculation using 
model NSC97f for the A7V interaction is exhibited in Fig. 
|l| as function of the number N of basis functions. The 
corresponding curve, marked 'Apn', gives the Bja value 
with respect to the horizontal straight line marked 'A-|-d 
threshold'. The 'Apn' asymptote serves then for defining 
the lowest particle-stability threshold, that of A-f ^H, in 



[ANN)s^i + A , [ANN)s^^ + A , (AAiV)^^ 1 -I- N 

(3) 

for 3-1-1 breakup clusters, and 



(AA)s=o + {NN) 



{AN)s + {KN). 



(4) 



with (S',S")=(0,l)-h(l,0) and (1,1) for 2+2 breakup 
clusters. We find invariably that the three rearrangement 
channels in which the two nucleons belong to the same 
d-like cluster dominate in actual calculations. This obser- 
vation, apparently, could justify the use of a AAd model 
for AA-H- However, as we shall see and discuss below, the 
results of using such a three-body model differ radically 
from those of the full four-body Faddeev- Yakubovsky cal- 
culations which retain the proton and neutron as dynam- 
ically independent entities. 



III. RESULTS 



Using the AA interaction which reproduces BAhipJi^^) 
(see Table 0) our calculations yield no bound state for 
the AApn system, as demonstrated in Fig. |^ by the lo- 
cation of the 'AApn' curve above the horizontal straight 
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FIG. 1; Convergence of Faddeev-Yakubovsky calculations for 
the binding energy of the Apn [S=l/2), AAd and AApn {S=l) 
systems with respect to the number A'' of basis functions. Val- 
ues of i?cutofF = 30 fm for AApn and Apn, and 60 fm for AAd 
were used. The A A interaction is due to Table 



FIG. 2: i3AA(AAH) calculated in a three-body AAd model as 
function of the scattering length aaa, for two exponential Ad 
potentials corresponding to versions / and /' in Table [II 



The solid squares correspond to the 'aaHc' AA interaction of 
Table The straight lines are drawn only to lead the eye. 



line marking the 'A + threshold' In fact our 

Faddeev-Yakubovsky calculations exhibit little sensitiv- 
ity to the strength of the AA interaction over a wide 
range, including much stronger AA interactions such as 
ND and ESCOO the latter one reproducing the (exces- 
sive) Baa value reported for the 'old' aaHc event For 
these AA interactions we get a bound only if the AiV 
interaction is made considerably stronger, by as much as 
40%. With four AiV pairwise interactions out of a total 
of six, the strength of the AiV interaction (here about 
half of that for NN binding the deuteron) plays a major 
role in the four-body AApn problem. In passing we re- 
mark that this is also apparent from the bounds derived 
in Ref. 
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for the four-body bound-state problem. Put 
differently, we know of no few-body theorem that would 
imply, for essentially attractive AA interactions and for a 
non static nuclear core d (made out of pn in the present 
case), the existence of a AAd bound state provided that 
Ad is bound. It is a remarkable outcome of the complete 
Faddeev-Yakubovsky scheme for four particles that such 
a natural expectation can be refuted by a specific calcula- 
tion. However, for a static nuclear core d, and disregard- 
ing inessential complications due to spin, a two-body Ad 
bound state does imply binding for the three-body AAd 
system ||2^. A discussion of the formal relationship be- 
tween these four-body and three-body models which do 
not share a common hamiltonian is deferred to a subse- 
quent publication. 

Our AAd model for aaH uses the AA interaction 
marked '^^He' in Table O plus Ad interactions that re- 
produce the low-energy parameters of the Apn Faddeev 



calculation specified in Tabic III. The dependence on 



the functional form chosen for the interpolating Ad inter- 
action potentials proved relatively mild. The results of 



such a AAd three-body Faddeev calculation using model 
NSC97f for the underlying AA^ interaction are shown in 
Fig. |l| as function of the number N of basis functions 
used in the expansion of the Faddeev components. The 
asymptote of the curve marked 'AAd' is now located be- 
low the horizontal straight line for the 'A+ threshold', 
so is particle stable. The figure may suggest that a A 
in is less bound, by about 0.1 MeV, than a A in 
( which in model NSC97f is bound by about 0.2 MeV). 
However, the {2J + 1) spin-averaged effective -Ba(^H) in 
is only Ba = 0.07 MeV and, since Baa 0.3 MeV, 
we have Baa > 2-Ba which is equivalent to stating loosely 
that the second A in aa^ bound even more strongly 
than the first one. This holds also for model NSC97e and 
it is a general property of the Faddeev calculation ^ . 

In Fig. 

within this AAd Faddeev model as function of Vaa (quan- 
tified by the value of the scattering length aaa) for two 
exponential Ad potentials corresponding to versions / 
and /' of model NSC97 (see Table |l|). The roughly lin- 
ear increase of i3AA holds generally in three-body AAC 
models (C standing for a cluster) over a wide range of 
values for Vaa §• For values Baa < 0.2 MeV, aa^ 
becomes unstable against emitting a A. This onset of 
particle stability for aa^ requires a minimum strength 
for the AA interaction which is satisfied for our choice 
of aaHs ^s ^ normalizing datum. It is also seen from 



we show Baa values calculated for ^^H 



the figure that the uncertainty in the location of aH(| ) 
bears serious consequences for the predicted binding of 
aaH; this is a particularly strong effect as the |^ state 
crosses the A -I- d threshold. Yet, we would like to em- 
phasize that no such sensitivity emerges within a genuine 
four-body model calculation which does not bind aa^ 
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long as the AN interaction is of the size constrained by 
single-A hypernuclear phenomenology. 

IV. DISCUSSION AND CONCLUSIONS 

In cluster models of the type A AC and AAC1C2 for 
heavier AA hypernuclei, where the nuclear-core cluster 
C=Ci+C2 is made out of subclusters Ci and C2, the 
ACj interaction (normally producing bound states) is 
considerably stronger than for AiV. Our experience with 
Faddeev-Yakubovsky calculations for xa^c Q, viewed 
as a four-body AAaa system, is that the relationship 
between the three-body and four-body models is then 
opposite to that found here for ^^H: the AAC1C2 calcu- 
lation under similar conditions provides higher binding 
than the AAC calculation yields. The mechanism be- 
hind it is the attraction induced by the AC1-AC2, AACi- 
C2, C1-AAC2 four-body rearrangement channels which 
include bound states that have no room for in the three- 
body AAC model. The binding energy calculated within 
the four-body model increases then 'normally' with Vaa- 



In conclusion, we have provided a first four-body 
Faddeev-Yakubovsky calculation for using NN and 
AA^ interaction potentials that fit the available data on 
the relevant subsystems, including the binding energy of 
^H. No bound state is obtained for over a wide range 
of A A interaction strengths, including that normalized to 
reproduce the binding energy of aaHc. We have traced 
the origin of this non binding as due to the relatively 
weak A7V interaction. This is in stark contrast to the 
results of a 'reasonable' three-body AAd Faddeev calcu- 
lation that binds j^^H provided the AA interaction is not 
too weak, say with — caa > 0.5 fm. 
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